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ABSTRACT 


Refinements to the bow and stern seal subroutines of the 
XR-3 Loads and Motion computer program are made. New lift 
fan performance maps are developed and presented. Comparison 
studies of simulated performance of the original XR-3 model 
with the old and new fan maps are made under varying 
sea conditions. Measured roll and pitch transient 
behavior of the XR-3 test craft in calm water is compared 
with the simulated response of both the original and the new 
model. In addition, simulated performances of the original 


and new models in irregular ahead seas are compared. 
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A. BACKGROUND 

. The speed limitation of the conventional displacement 
vessel has been recognized for many years. In an effort to 
effect a quantum increase in surface vessel speed the United 
otates Navy has been developing various craft whose principal 
means of support is other than hydrostatic lift. 

One type of craft currently receiving much attention is 
the surface effect ship or craft. In this category are the 
air cushion vehicles and captured air bubble craft (CAB). 
The general nature of these craft and their construction is 
Tew presenvedr by NObert lb aTrillocin Ret. ila In these 
craft either all or a major portion of its support is derived 
from a pressure differential between the atmosphere and a 
plenum chamber which iS open at the bottom. The speed 
advantages gained by this type of vessel are derived from 
two principal characteristics. It does not waste energy 
displacing a large volume of water and the frictional forces 
between the hull and water are greatly reduced by the craft 
riding on a cushion of air with only a minimum cf the craft 
structure actually in contact with the water. 

Vessels of this type are generally separated into two 
categories — "Air Cushion Vehicles" whose weight is entirely 
supported by the pressure differential in the plenum chamber 


and "Captured Air Bubble Craft" whose weight is partially 


diel 





| 
supporved by sidewalls which extend into the water. It is 
the second type of vessel in which the United States Navy 
has become primarily interested for ocean going applications. 

The term Captured Air Bubble (CAB) is possibly a enaner 
since in this type of craft the air does leak out and has to 
be continually replenished by supply fans. The leakage 
though is relatively slow when it is compared to the air 
cushion vehicle which has continuous air gap around the 
periphery of the craft and hence the fans that replenish the 
| air in the plenum chamber have to be larger and more powerful 
than those of a CAB of Similar size. 

This thesis is concerned with one aspect of the develop- 
ment of the Captured Air Bubble vehicle and that is the 
SVU atlOnwor such Crary Utilizing a dieital computer, 
specifically the Loads and Motion Program developed by 
Oceanices. ine. 

The basic rigid body analysis and spatial relationships 
of the Loads and Motion program are well documented in 
Ref. 2. Also covered are the principal dynamic and static 
approximations used in developing the equations of motion 
for the craft in its six degrees of freedom. 

A copy of this program was delivered to the Naval 
Postgraduate seneel tia Ocueper Of 19/2 whieh had been 
Specifically tailored for the Bell (100B) 100 ton craft. 

Leo and Boncal in Ref. 3 succeeded in converting this program 


to represent the XR-3 test craft. 





B. OBJECTIVE 

The purpose of this thesis was to take a second more 
detailed look into all of the aspects of seal modeling and 
the fan performance map with the idea of possibly making 
improvements in these representations before testing the 


computer model of the XR-3 under dynamic conditions. 
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II, FAN PERFORMANCE MAP ANALYSIS AND COMPARISON 


A. INTRODUCTION 

As the surface effect ship is predominately dependent 
Cmerne peenum lilving@eerce for its performance, the 
emaracveristics of the lifting fans are important if valid 
studies of the craft in question are to be made. 

Leo and Boneal reported in Ref. 3 that the steady state 
operating conditions of the XR-3 were obtained by increasing 
the leakage area of the craft until the proper equilibrium 
pressure was maintained in the main plenum. This area was 
twenty-six square feet. As the cross-sectional area of the 
plenum is less than twenty square feet this leakage area was 
out of proportion to properly represent the leakage area of 
the XR-3. 

While operating above hump speed the XR-3 test craft has 
been observed to have a gap underneath the sternseal of 
from four to five inches. Slight leakage does occur at 
the bowseal as reported by Ref. 3 but the total leakage area 
is estimated to be not over four square feet. 

To appreciate the effect this leakage area has on computed 
craft performance an understanding is required of the method 
by which the pressure is calculated in the plenum. 

When the program is executed an estimation of the initial 
operating conditions must be Supolredsby the user. Besides 
Sloe eolia thrust, initial conditzons for plenum pressure, 


displacement and pitch angle must be supplied. From the 
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displacement, pitch angle and velocity for calm water 
conditions, a plenum volume is computed. Using the leakage 
areas supplied, exit flow rates are computed for the bow 
Pies Ler) cealo. tne fan Subroutine is then called using 
the pressure in the plenum, bow and stern seals as primary 
arguments. Using these pressures supply flow rates are 
established for.the fans that Supply air to the plenun, 

and both seals. As the seals discharge into the main plenum 
Meme un Omeypnese Tagelow rates €stablish the total influx 
rate. The net flow being the difference between the influx 
and efflux rates establishes the new plenum pressure for. 
the next calculation interval based upon the adiabatic 
perfect gas law. This new pressure creates a new plenum 
lifting force which is then balanced by the defining differ- 
ential equations of motion resulting in a new operating 
condition of craft attitude and displacement. 

Meshes ChatwenmOtLOn in calm water the leakage area 
is of little or no consequence Since no additional leakage 
areas are normally developed beneath the sidewalls or seals, 
and the rate of change of plenum volume is relatively slow. 
On the other hand, when the program is executed in simulated 
sea states where gap areas are generated beneath the seals 
and sidewalls and plenum volume does change quite rapidly 
the leakage areas and fan maps become of prime importance. 

It must be emphasized though that it is not the leakage 
area itself that affects craft performance. It is rather 


the relative size of the fixed leakage area with respect to 








the variable leakage areas generated around the periphery 
of the craft by wave motion and craft attitude and position 


which is important. 


B. FAN MAP ANALYSIS AND DERIVATION 

Isolation of the factors which caused the large 
discrepancy between observed steady state gap area and 
calculated area required an examination of both the efflux 
and influx equations for possible errors. 

The equation used to establish leakage flow rates under 
the seals and the sidewalls are the same. This equation 
is given below 

Q = C pA, (2Ap/p) 1/4 

It is the general flow equation for a sharp edged orifice 
in turbulent flow conditions. As the orifice coefficient, 
Ces is the only variable in the equation that could be 
dependent on craft geometry and its range was limited between 
unity and six tenths; the efflux equations were eliminated 
as the source of the error. The only area left for possible 
error was the fan maps themselves. 

The fan maps in question are arrays which give fan flow 
rates in cubic feet per second as a function of plenum or 
seal pressure. They are entered as data in "INCON" as 


characteristics of the specific craft being modeled. 
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It was found that the fan maps constructed by Leo and 
Boncal to represent the fans installed in the XR-3 were Bees 
upon incorrect information. Actual fan performance data for 
the fans currently installed is presented'in Figure 1, and 
physical characteristics in Figure 2. 

The fan performance map given by the manufacturer, Joy 
Manufacturing Company, is based upon flow that is not impeded 
by air ducts. It is therefore optimistic for any system 
that contains long ducts as these ducts will modify 
considerably the fan performance by creating ducting losses. 

The best way to obtain the actual fan performance maps 
in a duct system would be to actually measure the flow rates 
as the pressure in the receiving reservoir (in this case the 
plenum).is varied. As this method was impractical in the 
case of the XR-3 another neteod of estimating the modified 
fan performance had to be devised. 

The method outlined below was utilized and although it 
probably does not give exact results the fan maps thus 
achieved should adequately represent the actual fan performance 


of the XR-3. 


lie fan data was provided by Associate Professor of 
Peremwory Do. Layton who is in charge of the actual XR=3 
test program at the Naval Postgraduate School. He has 
verified that the fans installed are Joy Model AVR87-55D2058 
Zacmmeemecie Gdava given in Figures 1 and 2 are correct. 
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FIGURE 1 
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PERFORMANCE AND FAN Map Curves FoR THE Joy AxIVANE FAN 
(MopeL AVR 87-55D2053) 
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The first step used in arriving at the new fan maps was 
to establish a realistic leakage area. Once this leakage 
area was sates a steady state efflux rate was 
established. Because in steady state operating conditions 
the influx and efflux rates must be equal the flow rate 
Supplied by the lift fans in total was established. It was 
at this point that a second assumption was made. This 
assumption was that there were greater losses in the ducting 
Supplying air to the bow seal and main plenum chamber than 
there were in the ducting supplying the rear seal. This 
assumption was based upon the fact that the ducting to the 
bowseal and main plenum had to carry the air forward 
approximately sixteen feet to their discharge orifices 
whereas the stern seal ducting was only a few feet long with 
orev one curve in ieeco impede fiow 

AS e@ result the fraction of total@iwow assigned Toemmec 
stern seal was greater on an individual fan basis than was 
the air supplied by the bow seal or main plenum fans. 
Another assumption was that the same map could adequately 
represent both the bow seal supply fans and those supplying 
the main chamber. This last assumption was based on similar 
Oieteeneems, andvconstruction. 

Now knowing the individual fan flow rates and the actual 
plenum pressure a pressure drop in the duct for steady state 
conditions was established. This pressure drop was the 
difference between the known plenum operating pressure 


(24.84 lbs/ft*) and the indicated fan discharge pressure 





for that flow rate. This pressure drop was then scaled 
for different flow rates based upon the fact that pressure 
losses vary with the square of the flow rates. In this 
manner new fan maps were computed for the bow seal and 
main plenum, and the stern sea. These maps are presented 
in Figures 3 and 4. 

The last step was to go back to the leakage area and 
make the slight adjustments needed to give the required 
plenum steady state pressure of 24.84 lbs/ft*. The final 
leakage areas established were .08 rts for the bow seal 
amd 3.79 ft* .for Uiomotemaesce tl, 

In the original fan subroutine the horsepower generated 
by the fans is computed based upon the infiux rate and the 
pressure differential between the atmosphere and the plenum. 
As the actual work done per unit time by the fans this 
Pipurescan be misinterpreted as the actual power needed 
by the fans to generate printed work rate. This would only 
hold true if the fans operated with one hundred percent 
efficiency and there were no ducting losses. In actuality 
fans of this nature which supply relatively large volumes 
of air at low pressure rarely operate at over seventy percent 
efficiency and can achieve that only over an extremely 
narrow operating range. These physical characteristics of 
fans were well presented by R.D. Moyer in 1946 in Ref. 4. 

For the fans that are used in the XR-3 the Joy Manufac- 


turing Company has provided the fan's horsepower requirements 
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on the same performance graph as the fan maps. This power 
requirement is plotted for the entire operating range of 

the fans and is fairly consistent at about 2 horsepower. 
This horsepower requirement does not change under the new 
fan maps since the fan map discharge pressures at the fan 
exhaust remains the same. This horsepower map was entered 
in the fan subroutine as a table of look up and actual power 
requirements were computed from it using flow rate as the 
entering argument. 

In the COLFIL subroutine (the routine that prints out 
all data and draws all graphs) the fan power label was 
changed to "fan power (ideal)" to avoid possible confusion. 
In the fan subroutine the print option was modified to give 
the actual power required as well as the ideal power. Also 
the fan efficiency is computed and is given as an output of 
mons Subeoutane. 

C. SIMULATED PERFORMANCE COMPARISONS BETWEEN 

ORIGINAL AND NEW FAN MAPS 

As stated earlier it was believed that changing the 
leakage areas and the fan maps would have a Significant 
effect only when the computer model was run under conditions 
where the plenum volume underwent comparatively rapid 
fluctuations and/or the leakage areas around the periphery 
Of eheenull varied. 

These assumptions were confirmed by Simulated runs under 
varying sea conditions. In all these comparison tests the 


two XR~3 models used were identical except that the original 
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fan maps were replaced by the newly developed ones in one 
model. Both simulated craft were run using the original 
seals as developed in Ref. 3. Craft weight was 5813 lbs, 
center of gravity 9.57 feet from the stern, speed was held 
constant at 20 knots. All runs were made into ahead seas. 
The sea conditions chosen were irregular sea state 1, 
irregular sea state 2, regular seas 1 and regular seas 2. 
The latter two were so named because the amplitudes and 
frequencies chosen for these regular seas runs were the major 
component of irregular sea state 1 and irregular sea state 2 
respectively. The actual amplitudes and frequencies for 
Piescurwnosarc vabulaveq on, page 27. 

In the simulation runs into regular seas 1 dramatic 
differences were noted in the fluctuation of plenum pressure 
and heave accelerations. With the original fan maps, 
variations in plenum pressure were limited to about 2.2 lbs/ft* 
whereas with the new maps the variations were over 6.9 lbs/ft°. 
This represents an increase in fluctuations of over 315 
percent. Because the plenum normally provides over seventy 
Mereent Cf the Supporting foree for the craft similarly 
large variations were noted in heave acceleration both at 
the bow and at the center of gravity. Pitch motion, though, 
did not change much being less than .1 degree in both cases. 

Results for the runs in other sea conditions are 
tabulated on page 28. The differences noted for the runs 


in irregular sea state 2 and regular seas 2 were not as great. 





This nas been attributed to the craft's ability to conform 
its attitude to the wave slope in the case ofthe longer 

Wave lengths associated with greater sea states. This ability 
to conform then reduces the high frequency volumetric changes 
in the plenum chamber, thereby reducing pressure fluctuations. 
Paesures 5 through l2 are the results of running both models 
in the simulated sea state one. By visually comparing the 
corresponding variances of the two models it is readily 
apparent that the fan maps do significantly change the 
performance characteristics of the model. By examining the 
table on page 28 it is also apparent that one must be careful 
in predicting behavior as the variation in measured variables 
does not vary directly with the sea conditions. The new fan 


maps have in effect defined a similar but different craft. 
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TABLE I. SEA STATE WAVE COMPONENTS 


Run 1 Regular Seas 1 Heading 180 


Wave Components Amplitude (ft) 


Frequency (rad/sec) 


1 0598 4.1526 


Run 2 Regular Seas 2 Heading 180 


Wave Components Amplitude (ft) Frequency (rad/sec) 


al wee. 1.8447 


Run 3 Irregular Sea State 1 Heading 180 


Wave Components Amplitude (ft) Frequency (rad/sec) 


Run 


y 


1 POOLS 2.553 

2 ose 3.256 

3 .0599 NAL 26 
i 70550 5.2961 
5 .O401 6.7545 
6 m6 2 8.6145 
fi .0165 10.9867 
8 0102 ae oes 


Irregular Sea State 2 


Wave Components 


Heading 180 


Amba. tude mort ) 


Frequency (rad/sec) 


1 .0168 1.1786 
2 L464 1.4745 
3 Nee 1.8447 
if -Achels) 2 SOT 
5 ne Ope Zoe talk 
6 5 vate Bj. eal iis 
7 0928 4.5185 
8 0600 Bob ee 
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TABLE II SIMULATED XR-3 PERFORMANCE WITH DIFFERENT 
FAN MAPS 


Comparison of XR-3 Loads and Motion Computer Program 
performance with original for maps with XR-3 with new 
fan maps are given below: 


RUN 


Regular 
seas 1 


Irregular 
Sea State l 


Regular 
Seas 2 


lrregular 
sea State 2 


QUANTITY 


Plenum Pressure 
Inakiycla, jtaoals 

Draft 

C.G. Acceleration 


Plenum Pressure 
Pitch Angle 

Draft 

C.G. Acceleration 


Plenum Pressure 
Pitch Angle 

Dice ae 

C.G. Acceleration 


Plenum Pressure 
Pitch Angle 

Draft 

C.G. Acceleration 
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FIGURE 5 
Computer SIMULATION PLENUM PRESSURE VERSES TINE, SPEED 20 KNoTS 
ORIGINAL FAN Maps, SEA STATE l 
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FIGURE 6 
Computer Simucation PitcH ANGLE VERSES TIME, SPEED 20 KNOTS 
ORIGINAL FAN Maps, SEA STATE 1 
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ORIGINAL Fan Maps, Sea State 1 
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FIGURE 8. 


Computer SIMULATION C.G. ACCELERATION VERSES TIME, SPEED 20 KNOTS 
ORIGINAL Fan Maps, Sea State 1 
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New Fan Maps, Sea State 1 
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III. SEAL MODEL REFINEMENTS 


A. INTRODUCTION AND BASIC MODELING TECHNIQUES 

Through tests conducted aboard the XR-3 by the NPS test 
group headed by Associated Professor D.M. Layton of the 
Aeronautics Department it has been demonstrated that the 
forces generated by the seals have an appreciable effect on 
tee overall perrormance Of the test crait. Sy varying the 
lower travel limit of the bow seal and contouring its shepe 
WeewslCady sstarc Patten aagle Gan be changed @s well as the 
thrust versus velocity characteristics. Although not as 
pronounced, changing the support cable length of the stern 
Se@ll BLO sil Wereetetsl toige nie jerealigl oe etslakjealvoyg sip 

Because of the demonstrated steady state pitch sensitivity 
of the XR-3 to seal position limits as set by the Support 
cables it was decided that more detailed modeling of the 
seals was needed before the overall response of the L&M 
computer model would reflect the behavior of the XR-3 under 
varying conditions of support cable length. 

The original XR-3 seal modeling as presented in Ref. 3 
was patterned after the rear seal in the L&E computer model 
of the Bell (100B) 100 ton CAB. This seal is of the "soft" 
pneumatic type and its representation is wejl covered by 
Kaplan, Bentson and Sargent in Ref. 2. 

The pneumatic "soft" seal model is based on the assumption 
that this type of seal is unstiffened and non-rigid. The 


mass of the seal was considered so insignificantly small as 


ou 





compared to the forces generated by the seal when it is in 
contact with the water as to have negligible effect on overall 
craft dynamics. The lifting forces generated by the seals 
were calculated by multiplying the wetted contact area of 
the seals by the pressure differential between the seal and 
the atmosphere for the bow seal and the plenum chamber and 
See in the case of the stern seal. 

Drag forces on the seals were considered to be due only 
to skin friction between the water and the seal. The equation 


used to compute these forces was 


™~ 


2 
eee 
ke 5 C 


pe gn 
ea di seal 


where p 1s the density of water, u the velocity of the craft 


and A. is the wetted surface of the seal. Car is the 


eal 
drag coefficient and is calculated using the Schultz-Grunow 


formula 


S = -427/ (log, oRe -- 407224 


GH: 


where Re is the Reynolds number computed based on the average 
wetted length of the seal and the Reynolds number is 


calculated by: 
u x (wetted length) 


ne <6 
ree 10 


where LAS aie ft“/sec is tne kinematic viscosity olf the 


water, 
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The actual construction of the seals is well described 
tener, 3. Each weneconsilete of two pneumatic lobes 
separated by a common membrane. The membrane has limber 
holes in it to equalize the pressure between the lobes. The 
lower surface of the bottom lobe is shaped by twelve steel 
Baevens that aet as stiffeners. These battens are spaced 
uniformly across the bottom surface at right angles to the 
hinge line. It is this surface that comes in contact with 
the water, and will hereafter be referred to as the "leading 
face’ of the seal. 

The forward edge of the leading faces is hinged near .the 
top of the plenum chamber allowing the seal to pivot about 
this line when waves are encountered. 

The pressure in the seals iS maintained by axial fans. 
The alr is delivered by ducts that discharge into the top 
er the upwers lobe. The seals in turn vent to the main 
plenum by means of ducting. Back pressure in the vent 
ducting can be controlled by butterfly valves that can be 
manually adjusted. Normally the seals are operated so that 
iney Matitain about F=2 PSF more pressure than the plenum. 

Downward travel limits of the seals are controlled by 
adjustable support cables attached to the rear edge of each 
steel batten at the trailing edge of the seal's leading face. 
In addition a second support cable is attached to each 
stiffener apvoroximately midway between the hinge line and 
the trailing edge. Figure 13 shows a cutaway of the seals, 


only two cables are drawn in for clarity. 
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Both the bow and the stern seal extend the width of 
the plenum chamber and are identical in construction. The 
leading faces are approximately 49 inches long from hinge - 
line to trailing edge. The perpendicular distance from 
Sidewall keels to the hinge line is 22.5 inches. 

As stated before the original modeling of the Seals did 
Meme eanc inlbo account any of the seal's ability to flex or 
peevot Ol its hinge. In addition no Aevanot was made to 
model the support cables and their ability to affect craft 
attitude. In the new modeling the support cables have been 
icorporated. In the stern seal modeling only the rear 
Support cable is modeled since this seal rarely comes in 
contact with the water. In the bow seal subroutine both the 
rear support cables and the center support cables have been 
included since this seal does contact the water. Before any 
attempt at reviSing the seal Subroutines could be undertaken 
it was necessary first to study the action of the seals while 
Peemeratt Was actually underway. 

The semi-rigid construction of the seal dictated a highly 
complex distributed parameter model of the seals be developed 
if the seals were to be fully modeled without any approxima- 
puemo we tn addition to the fact that practically all of the 
physical characteristic parameters of the seal such as spring 
Somenemia tor the steel stiffeners and fabric flexibility 
were unknown such a model would require an inordinate amount 
of calculation time in the computer. Since extra computation 


time was to be avoided simplifying assumptions had to be made. 
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To be accurate these assumptions or approximations had to 
Deebeeed On actual observations of the seal or on previously 
made approximations whose effect has already been investigated. 

In Ref. 2 the approximation of eliminating seal dynamics 
and ae setting of seal pressure at a constant level above 
plenum pressure has been thoroughly investigated and was 
incorporated in the original model of soft pneumatic seals. 
These approximations have been continued in the present seal 

modeling. 

| The other modeling approximations made were based upon 
actual observations made of the seals while the XR-3 was in 
mMeGLOn. These observations are described below: 

1. For small water perturbations (approximately 3 inches) 
the bew seal flexed and appeared to pivot slightly about the 
center support cable causing the angle at the hinge line 
Ie slightly increase, but almost insignificantly. 

2. For longer perturbations the seal both flexed as 
described above and then began to pivot on the hinge. After 
the seal rose about six inches no additional flexing took 
place. 

3. The bow seal maintained contact with the water at all 
times and never gaped because of an inability to follow wave 
moO waEON . 

4. The bow seals penetration of the water was practically 
independent of water motion in relation to the craft, but 


appeared to be inversely related to the craft's velocity. 
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When the XR-3 was at rest penetration was about five to six 
inches and when traveling at 20 knots appeared to be 
approximately 1 inch. 

5. The stern seal rode out of the water while the craft 
peerated above "hump" speed. The height of the air gap thus 
created appeared to remain almost constant at four to five 
inches. Water rarely contacted the seal. 

Using the above described observations and the preceding 
Simplifying approximations concerning seal dynamics and 
pressures both the bow and stern seal were modified to 
meamlece tae actaicn Of the actual seals. 

The approximation given in Ref. 2 of equating the wetted 
‘length of the seal with the water height above the sidewall 
Weel at ths Sseal's trailing edge location was eliminated. 

In its stead a two-dimensional "table of look up" was 
developed which gave wetted length as a function of the 
immersion depth of the seal's trailing edge in the water and 
the seal's trailing edge deflection height above the keel at 
the trailing edge location. In other words if the seal was 
deflected upwards because of wave motion or pitching motion 
of the craft the deflection height of the trailing edge 

would increase and the seal's leading face would become more 
parallel to the water's surface thereby increasing the wetted 


area of the seal. 





The "table of lookup" was constructed using the 
information obtained from a graphical deflection model 
similar to the one presented in Figure 14. In this model 
the leading face of the seal is represented in varying 
amounts of deflection height from full down (point A) to 
full up against the seal framework (point B). Full down 
is zero deflection anCmeomumecre Le = tira mememcane Of the 
seal is even with the bottom of the sidewall keel, full: up 
is maximum deflection. This deflection height was divided 
imo Sesments Of One winch for ease of computation. Points 
veGwe Marked upon each leadime face line which represented 
the surface of the water on the leading face under varying 
depths of trailing edge immersion. The distance along the 
meading face to the trailing edge from these points then 
represented the length of the leading facé that was immersed 
in the water or the "wetted length". The points for similar 
immersion depths were connected for the varying seal deflec- 
Peon ftetaiuemanc these construction lines are labeled with 
their corresponding immersion depths in inches along the 
zero deflection height leading face. The construction line 
corresponding to an immersion depth of zero inches is, of 
course, congruent with the trailing edge of the seal since 
if the seal did not penetrate the water to at least some 
degree there could be no "wetted length" unless the leading 
face was exactly parallel to the water's surface, which 
never occurs. These construction lines were drawn for 


immersion depths of the trailing edge of 0 to 5 inches. 
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The actual immersion depth of the trailing edge is 
calculated based upon the observation that at rest the 
immersion depth was approximately 5.5 inches and while 
traveling at 20 knots was estimated at a little less than 
one inch. Since the hydrodynamic force on the seal is 
proportional to the square of the craft's velocity the 
immersion depth was assumed to be approximately inversely 
proportional to some factor of the square of the speed. 

Using 5.5 inches as the known immersion depth at zero 
meeced and 2.75 inehes at cO knots, the following empirical 


formula was derived which gives reasonable results 
Demme 555/Get =) 
' ° 25) 


where D(in.) is the depth in inches and u is the velocity 
in feet per second. 

Using the calculated immersion depth and the’ seals 
calculated deflection height the table is entered to give 
the wetted length. This table of look up is presented 
graphically in Figure 15. The non-linear nature of the 
Cilcvesmis readily apparent. An important aspects Ole this 
tbaeilewapparent in Figure 15 1s that within abour the first 
ten inches of seal deflection the wetted length of the seal 
remains fairly constant, especially if the craft is operating 
at speeds where the immersion depth is approximately one inch. 
This figure also shows that the maximum wetted length is 


equal to the length of the leading face (just over 48 inches). 
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As the seal deflects and becomes more horizontal the wetted 
Surface area becomes greater given a constant immersion 
Gepc h. 

This same "table of look up" is provided in the stern 
seal modeling but because the seal does not ride in the water 
an assumed "depth" had to be used as an entering argument. 

ae assumption was made that since the stern seal does 
not ride on the water its lifting force has to be the result 
of some form of air pressure. This force could only be 
generated by the impingement of air particles on the sSeal's 
feading face as they rushed to escape underneath the seal. 

At this point an assumption was made that the effective 
wetted area occurred where the kinetic energy in the air mass 
was approximately one quarter that of the air escaping under- 
MeawiemunMemoacal., 2s the air in the plenum is at such a 
pressure that it can be considered incompressible this would 
Occur where the cross-sectional area was approximately two 
times the cross-sectional area beneath the seal. As the 
width of the plenum chamber is constant and the gap height 
is approximately 4.5 inches at all times the point on the 
seal where the cross-sectional area is two times the gap 
area is the point along the leading face where the height 
from the seal oe the water is two times the gap height or 
about 9 inches. Subtracting the gap height gives an 
effective "phantom" immersion depth of 4.5 inches; this was 


rounded to five inches. 
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Lift forces generated by the stern’seal were considered 
to be the product of the seal area found by using the above 
Phantom immersion depth and the seal deflection height, and 
the pressure differential between the stern Seal and plenum 
chamber. 

The above lifting force calculation has only minimal 
effect on simulated performance of the XR-3 Since the 
pressure between the seal and the plenum chamber is normally 
only one to two PSF. The drag forces on the stern seal are 
calculated in the same manner as they are in the bow seal 
Subroutine but are subsequently set to zero since the seal 


rarely contacts the water. 


B. BOW SEAL MODELING 

Contained herein is a detailed analysis of the actual 
bow seal subroutine. It is coded in standard FORTRAN. By 
reviewing the actual Subroutine the actual approximations 
used in the routine may be understood. A full listing of 
the BOWSL subroutine is contained in Appendix A along with 


the required entries in the data input statements. 


GAP(J) 


I) eh NT 
= ( 

) EL SKUC) = 

) 

) 

) 


Sue 
0.0 
WETLEN (J) OO) 
ELSKID(J) ore0 


1 CONTINUE 


(1 
(2 
(3 
(4 
cS 
(6 
The above statements (1 through 6) zero all variable arrays 


that will be used every time BOWSL is entered. GAP(J) is 


the height of the air gap at each station along the bow seal. 
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ELSKI(J) is the water height above the keel at the trailing 
edge of the seal. WETLEN(J) is the wetted surface length 
at each station. ELSKID(J) is the deflection height of the 


seal trailing edge measured from the botton of the keel of 


the craft. 
(7) ALBS = 0.0 
Coy Fx =o. 0 
CODY 6206 
Gley (FZ =". 0 
Gage = "0.0 
Gite = 0-6 
(13) FN = 0.0 


ine preceding seven statements initialize the single valued 
variables that are computed each time BOWSL is entered. 
ALBS is leakage area of the bow seal. FX is the forces 
parallel to the craft's axis of symmetry (centerline). 
FY is the horizontal forces perpendicular to the FX forces. 
mo 1s tne vertical forces™ FK, FM, FN are the moments about 
mee eaters! GxiS and X axis respectively. 

(14) DELPBG = PBS - PB 

(15) IF (DELPBG.LT.0.0)DELPBG=0.0 

Equation 14 computes the pressure differential between 

the plenum chamber (PB) and the seal (PBS). This differential 
is actually fixed in INCON as an entering argument. Equation 
15 assures that the seal differential pressure is never 
negative. This is reasonable since the seals vent to the 


chamber. 





(16) ARGO = ELMAXB/CORLEN 
Cr) ANGO = ARSIN(ARGO) 
(16) Xl = XBS+ZBS*THETA-CORLEN *COS(ANGO) 


The three equations above (16-18) establish the moment 
tiie ne = Circe enmaronme ne Craiiecentervof gravity 
to the trailing edge of the bow seal. ELMAXB is the rear 
cable length. CORLEN is the chord length of the bow seal. 
XBS is the distance from the center of gravity to the bow. 
ZBS is the vertical distance of the bow hinge line below 
the center of gravity. Theta is the pitch angle of the 


eran. 


(19) Zl = —~Z-ABS+XBS*THETA-ELMAXB*COS (THETA ) 
Equation (19) establishes the moment arm in the -Z 
Girection to the water's surface through which the drag 
forces act to create a moment about the Y axis. 
(20) DPTHFT = (5.5/(1.+(U/25.))¥**2.0/12.0 
Equation (20) calculates the immersion depth of the 
Seal Dasccmupon craft velocity, u. 
(21) IF(CENCAB.GT.1.1875) CENCAB=1.1875 


Logical statement (21) sets the maximum length of the 
center support cable to 3 inches more than the maximum rear 


GaekesVveneeh, ELMAXB. 
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(22) DO 3 K=1,N 


3 CONTINUE 


The above loop calculates for each seal partition point 
the water height above the keel ELSK1(K), the wetted length, 
WETLEN(K), and the seal deflection height ELSKID(K) in 


preparation to calculating the forces and moments. 


(23) DPFT(K) = DPTHFT 
Equation (23) initializes the immersion depth of the 
seal. 


(24) ELSK1(K)=(ETA(3,K)-DTABX(K) *(XX(3,K)-X1)-21) 
+YY(3,K) *PHI+XLXPWV#WATSLP 


Equation (24) establishes the water height above (or 
below) the keel at the bow. ETA(3,K) is the height of the 
water at the center of gravity above mean. DTABX(K) is 
the correction factor for moving this height away from the 
C.G. (Center of Gravity). (XX(3,K)-X1) is the distance the 
correction factor has to act on. 41 is the height of the 
keel in relation to the C.G. YY(3,K) is the displacement 
of the partition from the craft center line. PHI is the 
roll angle. XLXPNV*WATSLP is the correction factor for the 
water slope caused by the plenum pressure acting for a 


finite period. 
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(25) IF(ELSKI(K).GT.HINGHT) ELSKI(K)=HINGHT 


Equation (25) limits the water height to the hinge 
height of the bow seal. 
(26) IF ( (HINGHT-ELSKI(K)+DPFT(K)).GE.ELMAXB) 
DPFT(K)=ELMAXB-HINGHT+ELSKI(K) 
Equation (26) adjusts the immersion depth of the seal 


so that the seal can go no lower than that permitted by the 


rear cable ELMAXB. 


(27) Poin. e) DPFT=0.0 


DPFT is used in computing DPIN which is an entering 
argument into the wet length table BWSL. Therefore 
equation (27) limits this argument to a minimum value of 


Ber. 


(28) ELSKID(K )=(ELSKI(K)-DPFT(K))*12.0 
Equation (28) calculates the deflection height of the 
seal (ELSKID(K)) in inches based upon water height, ELSKI(K) 
and the immersion depth DPFT(K). 


(29) IF( (HINGHT-ELSKID(K)/12.0).GE.ELMAXB) 
ELSKID(K)=(HINGHT-~ELMAXB)*12.0 


Equation (29) limits the downward travel of the bow seal 
to ELMAXB the rear cable length and adjusts the deflection 


height ELSKID(K) to reflect this lower limit. 





' DPIN=DPFT(K)#12.0 


(30) 

(31) MM=DPIN 

(32) MM1=MM+t1 
(339) MM2=MM1+1 
(34) DINC=DPIN-MM 


The equations above, (30) through (34), establish 
parameters for entering the bow seal wetted length array, 
BOWSL. DPIN is the immersion depth in inches, MMl1 is the 
lower limit table calling integer, MM2 is the upper table 


calling integer and DINC is the incremental difference. 


(35) GAP(K) =-ELSKI (K)+ (HINGHT—ELMA XB) 


Equation (35) computes the air gap height based on 
water height ELSKI(K) which is positive upwards and the 


limit of the seal travel (HINGHT-ELMAXB). 


(36) ietcee,LT.O.0) GAP=0.,0 


Equation (36) limits the gap to positive values. 


(3a DAR EWeKTD(K).GE.0.0) GO™@O 2 
(38) WETLEN (K)=ELSKI(K) 
(39) GO TO 3 


Equations (37) through (39) set the wet length equal to 
the water height if it is negative. This computation is 
used in case one seal partition point has a negative value 


while the others on either side are positive. 





(40) MM3=ELSKID(K) 

(42) MM4=MM3+1 

(42) MM5=MM4+1 

(43) DLING=ELSKID(K)—-MM3 


Equations (40) through (43) establish the final two 
entering arguments for retrieving the four points needed 
from BOWSL to compute the wetted length. ELSKID(K) is the 
deflection height of the seal. MM4 is the lower entering 
integer based upon the deflection height, MM5 is the upper 


integer. DLINC is the incremental difference used for linear 


interpolation. 

(44) BWSL1=BWSL (MM1 ,MM4 ) 

(45) BWSL2=BWSL(MM1 ,MM5) 

(46) BWSL3=BWSL(MM2 ,MM4) 

(47) BWSL4=BWSL(MM2 ,MM5 ) 

(48) BWSLA1=(BWSL2-BWSL1) *DLINC+BWSL1 

(49) BWSLA2=(BWSL4-BWSL3) *DLINC+BWSL3 

(50) WETLEN (K )=( (BWSLA2-BWSLA1 ) *DINC+BWSLA1)/12.0 


Equations (44) through (50) identify the four points of 
the BWSL array which bracket the actual value of wetted 
Temeth. linear interpolation is used twice to find the 
average based on deflection height and once to find the 
final wetted length. The array BWSL was calculated in inches ; 
thus, to convert to feet the final average is divided by 


1250; 
cr) DO 8 J=1,N 


8 CONTINUE 
times Loop makes all the force calculations based upon 


wetted length and seal position calculated previously. In 


S18, 





addition it calculates the leakage area based upon the steady 


State leakage area and the variable gap heights. 


(52) WETLAV= (WETLEN (J+1)+WETLEN(J))/2.0 
(53) IF(WETLAV.LE..001)GO TO 6 

(54) DPFTAV=(DPFT(J+1)+DPFT(J))/2.0 
(55) ELSKIA=(ELSKI(J+1)+ELSKI(J))/2.0 
(56) ELSKDA=(ELSKID(J+1)+ELSKID(J))/2.0 


Equations (52) and (54) through (56) average the 
wetted length, immersion depth, water height and deflection 
height respectively, for each seal segment based on the 
Metues calculated for che stations on either side. 
Equation (53) routes the program around the remaining force 
calculations if the wetted length of the seal is zero or 


insignificant for calculation. 


Cay) SEALHT=HINGHT-ELSKDA 


Equation (57) is the calculation of the distance between 
the top of the seal's supporting frame and the bottom of 
the trailing edge of the seal. 


(58) . DIFF=2.0*CENCAB-(SEALHT+0.5) 
(59) IF(DIFF.GT.0.5) DIFF=0.5 


The DIFF equation (58) calculates the spring distance 
of the steel battens based upon the center support cable 
emer CeNCAB ana the position of the seal, SEALHT. 
BemieaGion (59) limits the spring movement to .5 feet or 
six inches, the point where rigid body analysis comes into 


full effect. 





(60) ARM1B(J )=X1+WETLAV/2.0 
(61) ARM2B(J )=ZS-ELSKIA+DPFTAV/2.0 


Equations (60) and (61) calculate the moment arms in 
the X and Z direction respectively based upon the length to 
the trailing edge of the bow seal, Xl, and half the wetted 
length, WETLAV/2.0, for ARM1B. ARM2B is based upon the 
height of the center of gravity ZS, the water height, 
ELSKIA and half the immersion depth of the seal, DPFTAV/2.0. 
Th wigieveGie©..25 ) GO TO 4 


2) 
3) DFBS(J)=-DELP*DELYBS*WETLAV 
4) GO TO 5 


ON ON ON 


( 
( 
( 


Equations (62) through (64) calculate the vertical force 
acting on the segment if the deflection difference is less 
twa .co meet Or three inches. If the deflection difference 
iomercaverevman .25 Or 3 inches the vertical force equavien 
has a spring constant term in it acting upon the moment 


summations about the hinge line. 


4 FORLEN=XBF-WETLAV 
IF (FORLEN.EQ.00) GO TO 5 
ARGW= (HINGHT-ELSKIA) /FORLEN 
IF(ARGW.GT.1.0)ARGW=1.0 
ANGW=ARSIN( ARGW) 
FORCOS=COS (ANGW) 
GO TO 6 

5 FORCOS=0.0 


LN LON LON LON LON ONO ON 
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Equations (65) through (72) calculate the length of the 
Seal not in contact with the water, FORLEN, and the angle 
which this length makes with the water ANGW. From the angle 
the cosine of ANGW is calculated, FORCOS, for future use 


Gaecalculating the vertical moments about the hinge line. 








(73)  DFBS(J)=-DELP*DELYBS*WETLAV-DELP*FORLEN *DELYBS 
* FORCOS * (FORLEN*0.5*FORCOS ) / 
FORLEN*FORCOS+WETLAV/2.0))# 
((DIFF-0.25)*4.0) 

The above equation calculates the vertical force acting 
upon the bow seal by summing vertical moments about the 
hinge line and equating it to the moment created by one 
vertical force vector acting on the midpoint of the wetted 
surface. The term "(DIFF-0.25)*4.0" introduces a spring 
Gistance of three inches through which the vertical force 
transitions from being only based upon the product of the 
wetted area and the pressure dagifferential to one based upon 
equalizing moments about the hinge line. After the seal 
bends this three inches it is asSumed to be completely rigid 
and therefore if the seal is not undergoing accelerations, 
the moments must balance about the hinge line. 
ARG=0.5*RHO*U*U *WETLAV*DELYBS 
RESKI=U*WETLAV/ENU 


ae 
a) 
6)  CDTSKI=0.427/(ALOG10(RESKI)-0.407)**2.64 
7) 


( 
} 
( TSKIB(J)=-ARG®CDTSKI 


I~ AI 


The above four equations calculate the drag forces of 
the seal based upon the previously discussed Schultz-Grunow 
formula. RESKI is the Reynolds number,CDTSKI is the drag 


coefficient and TSKIB(J) is the drag force on the segment. 


U1 
CO 





7 CONTINUE 


(78) FX=FX+TSK1B(d ) 

(79) FZ=FZ+DFBS (J ) 

(80) FK=FK+DFBS (J ) *YAVGB(J) 

(81) FM=FM-DFBS (J) *ARM1B(J)+TSKB(J) *ARM2B(J) 
(82) FN=FN-TSK1B(J) *YAVGB(J) 

(83) ALBS=ALBS+GAP(J)+GAP(J+1)) *DELYBS/2.0 


Equations (78) through (83) calculate the forces and 
moments for each segment and sums them to give total forces 
and moments. YAVGB(J) is the average distance off centerline 
of the individual segment. 

8 CONTINUE 
ALBS=ALBS+BLEAK 


(84) 
(85) 
(86) SQFAC=SQRT(2.#*ABS(PBAR)/RHOINF) 
ey) QL=CFBS*ALBS*SQFAC*SIGN(1. ,PBAR) 


The last equations calculate the flow rate of air based 
upon the leakage area ALBS and the pressure differential 
between the plenum chamber and the atmosphere. BLEAK is 
the seal leakage entered in INCON to account for slight 
leakage where the seals do not quite maintain contact with 
the sidewall. 

The remaining statements of subroutine BOWSL give the 
following values as printed output if the bow seal print 


switch is on, GAP, WETLEN, FX, FY, FZ, FK, FM, and FN. 


C. STERN SEAL MODELING 

The stern seal subroutine iS covered in the Same manner 
as the bow seal subroutine. The complete listing of the 
stern seal subroutine is contained in Appendix B along with 
iievrierions On the proper entries to be used on the data 


cards. 
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(1) DO 1 J=1,11 
(2) GAP(J)=0.0 
(3) ELSKI(J)=0.0 
(4) Eons )=0.0 
(5) AIRLEN(J)=0.0 
(6) 1 CONTINUE 


The above loop and equations zero all variable arrays 
that will change each time subroutine STNSL is entered as 


was the case in the bow seal subroutine. 


(7) EFDEP = 6.0 
Equation (7) establishes the effective depth of the seal 
at five inches. This figure is used as an entering argument 


in the stern seal contour mapping array CTNSL. 


(8) NESS) 
(9) eae 4 (0 
(10) IWR) 10 
(ia) FK=0.0 
(iz) FM=0.0 
(13) FN=0.0 
(14) AGAPI=0.0 
(15) AGAP2=0.0 
(key, AGAPA1=0.0 


The nine equations above zero all single valued variables 
that are computed each time STNSL is entered. 


(17) DELP=PSS-PB 
(18) IF(DELP.LT.0.0) DELP=0.0 


The two statements above fix the pressure differential 
between the plenum (PB) and the stern seal (PSS). Equation 
(17) is first used to arrive at the difference between 
these two quantities and the logical IF statement is used 


to insure all values are non-negative. As in the bow seal 
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subroutine this is logical since the air supplied to the 


stern seal vents into the plenum chamber. 


(19) PBAR=PB-PINF 


This equation establishes the pressure differential 
between the plenum ¢PB) and the atmosphere. It is used in 
computing leakage flow rates. 


C210) SINDIF=S INTH-COSTH*THETA 
C2b) COSDIF=COSTH+SINTH * THETA 


The above equations (20) and (21) establish the angle 


of the XR-3 with respect to the horizontal and the vertical, 


respectively. 
(22) X1=XSS+ZSS *THETA-XLF*SINDIF 
(23) Z1=(-Z-ZSS+XSS *THETA-ELMAXS*COS (PHETA) ) 


The above two equations establish lengths in the X and Z 
Girections which are used to calculate the water height above 
the keel at the stern seal location. 

(24) N=NSTA(4) 


(25) DO 2 K=1,N 
2 CONTINUE 


The above loop is used to establish the gap at each 
station and the phantom "wetted length" of the seal. 
NSTA(4) is the number of stations chosen for the seal and 


can be up to eleven, thereby defining ten seal segments. 


oak 





(26) ELSKI(K)=(ETA(4,K)-DETADX(K)*(XX(4,K)-X1) 
~Z1)+YY(4 jK) #PHI 
(27) IF(ELSKI(K).GT.HINGHT) ELSKI(K)=HINGHT 


Equation (26) computes the height of the water above 
the keel at the stern seal location. ETA(4,K) is a length 
correction factor that is dependent on seal shape. Since 
in the case of the XR-3 the seal is flat across the hinge 
line this term is zero. DETADX(K) is a wave correction 
factor that is multiplied by the seal distance from the 
center of gravity (XX(4,K)-X1) to give the wave height at 
the stern seal location. The logical IF statement limits 


the water height to the top of the seal. 


(28) ELSKIL(K)=ELSKI(K)+GPS 
(29) IF(ELSKIL(K).GT.HINGHT) ELSKIL(K)=HINGHT 
(30) IF(ELSKIL(K).LT. (HINGHT-ELMAXS) ) 


ELSKIL(K)=HINGHT-ELMAXS 


ELSKIL(K) is the height of the trailing edge of the 
Stern seal above the keel. It iS measured in feet. GPS alts: 
the air gap height that was observed to be constant under 
normal running conditions. It is .4644 feet. The logical 
IF statement (29) establishes the upper limit of _seal 
deflection where it is stopped by the seal Supporting 
framework. The second logical IF statement establishes the 


lower limit of seal deflection height. 


(31) GAP(K)=—-ELSKI(K)+(HINGHT-ELMAXS ) 
(32) eeecre Ghee. 0.0)GAP(K)=0.0 


The equation (31) and logical IF statement (32) establish 
the gap height below the lower stop limit set by the support 


cable length ELMAXS. 
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MM1=ELSKIL(K) *12.0 

MM2=MM1 

MM3=MM2+1 

DLINC=ELSKIL(K) #12. 0-MM1 

STNSL1=CTNSL (MM ,MM2) 

STNSL2=CTNSL (MM ,MM3) 
AIRLEN(K)=((STNSL2-STNSL1) *DLINC+STNSL1)/12.0 


i a aT an an al, Te alli 
Ud UO UO OO OO OD 
\O O~) OV Ee 
Ne Na a 


Equations (33) through (35) establish the entering 
arguments of the seal contour array CTNSL in the same manner 
as was done in the bow seal subroutine. AIRLEN(K) is the 
phantom "wetted length" and is given the name AIRLEN 
because the seal is supported by air. 

(40) N=NSTA(4)-=1 
(41) DO 5 J=1,N 
5 CONTINUE | 

The above loop accomplishes the actual force and 
moments of each seal segment and totals the forces and 
moments. It also computes the average gap height based upon 


the gap heights of the individual stations. 


(42) ELSKIA=(ELSKI(J+1)+ELSKI(J))/2.0 
(43) ELSKLA=(ELSKIL(J+1)+ELSKIL(J))/2.0 
Chip AIRLAV=(AIRLEN(J+1)+AIREN(J))/2.0 
The three above equations average the water height above 


the keel, ELSKIA;sthe average seal deflection height ELSKLA, 


and the average phantom wetted length ALRLAV. 


(45) AGAP=ELSKLA-ELSKIA 
(46) AGAPI1=AGAP 


Equation (45) establishes the average height of the air 
gap at the seal segment. AGAP1 is a variable used in further 


eabeulations. 
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(47) IF(AGAP.LT.GPS) AGAP=GPS 
(48) IF(AGAP1.GT.GPS)AGAP1=GPS 


The above logical IF statements establish the lower 
limit for AGAP at .4644 feet at the maximum value of AGAP1 
eee Ou feet. 


(49) ARM1S=XX(4 ,J)+ELSKIA/2 
(50) ARM2S=ZS-ELSK1A 


The above equations (49) and (50) calculate the moment 
arms through which the vertical lifting force acts to create 
the major component of the FM moment about the Y axis; the 
second equation creates the moment arm with respect to the 


face through which the drag forces act. 


(51) DFSS (J )=~DELP*DELYSS#AIRLAV/(GPS/AGAP) **2.0 


Equation (51) calculates the vertical force acting upon 
the seal and with the term (GPS/AGAP)**2,.0 allows for this 
force to decrease with the Square of the ratio of the actual 
gap width AGAP and the nominal width GPS. This is based on 
the assumption that the air exiting beneath the stern seal 
Supports the stern seal and that the liit forces generated 
by the exiting air are inversely proportional to their 


VieEMOCiLL YY. 
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(52) ARG=.5RHO*U*U* AIRLAV*DELYSS 

(53) RESKI=U* AI RLAV/ENU 

(54) GBESKI=.4277 GALOGIO(RESKI ) =. 407 )*##2064 
(55) TSKiS Gd )=—-ARG*CDTSKI 

(56) Skis (J )=000 


Equations (52) through (55) calculate the theoretical 
drag forces based upon the average phantom wetted length 
AIRLAV as if the seal was actually in contact with the water. 
fie ochultz Grunow formula is used as in the bow Seal routine 
to calculate the drag forces based upon the Reynold's number 
which is based on AIRLAV the average "wetted" length. 
Equation (56) removes this drag calculation since the seal 


is presumed to ride on a cushion of air and is hence 


negligible. 
(57) FX=PX+TSKIS (J) 
(58) PZ=FZ+DFSS (Jd) 
(59) FK=FK+DFSS (J) *YAVGS(J) 
(60) FM=FM-DFSS (J) #ARM1S(J)+TSKIS(J) *ARM2S (J) 
(61) FN=FN-TSKIS (J) #YAVGS(J) 


The five equations above calculate the forces and 
moments of each segment and sums them for all the segments 
to give the total seal forces and moments. 

ALSS=ALSS+ (GAP(J )+GAP(J+1) ) *DELYSS/2 .0 


(62) 
(63) AGAP2=AGAP2+AGAP1 
(64) AGAPA1=AGAP2/J 


The above equation (62) calculates the total gap area 
ALSS that is due to the seal reaching the downward limit 
of its travel and the gap thus increasing past its nominal 


value of .4644 feet. - Equations (63) and (64) calculate the 
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average gap height including the nominal .4644 height ana 
allows the nominal gap to close if the seal is slammed 


against its upper limit, the upper seal support framework. 


(65) ALSS=ALSS+ALEAK* (AGAPA1/GPS) 


Equation (65) calculates the leakage area under the seal. 
ALSS is the leakage area due to the gap width increasing 
beyond the nominal .4644 feet. The term ALEAK is the 
nominal leakage area of 3.79 square feet normally observed. 
Multiplying ALEAK by AGAPA1/GPS allows this nominal leakage 
area to decrease if the seal is limited in its upward travel. 


(Sia) SQFAC=SQRT(2. *ABS(PBAR/RHOINF) 
(67) QL=CFSS*ALSS*SQFAC*SIGN(1. ,PBAR) 


Equations (66) and (67) calculate the flow rate based 
upon the leakage area and the pressure differential between 
the plenum and the atmosphere. The equation assumes orifice 
type leakage. 

The remaining statements give print outputs of the terms 
Gee, ALRGEW, FX, FY, Pa, FK, FM and BN if the stern seal 


prim SWitch is en. 
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IV. DATA ANALYSIS OF ACTUAL CRAFT TRANSIENT 
MOTION AND COMPUTER SIMULATION 
A. INTRODUCTION 

In order to verify the modeling changes made in the 
Seals, data was collected on the XR-3 test craft for 
conditions of induced transtent motion in pitch and roll. 
To induce roll action the test craft was sallied by 
shifting body weight laterally near the C.G. until maximum 
Sscrtiiabions buiit Wp, baem the forcing action was stopped. 
For pitch oscillations the chase boat was veered across 
the path of the XR-3 to produce a wave front with its wake. 
By having the XR-3 pass over the wake it was intended to 
induce a transitory pitch oscillation. 

These attempts to produce a2 transient response cannot 
be viewed as entirely successful. The principal reason for 
the limited success was that it was extremely difficult to 
induce a large enough perturbation in either pitch or roll 
so that the normal background oscillations could be ignored. 
Figures 16, 17 and 18 are graphs of the roll and pitch 
or ons< of the XR-3 when making its twenty knot calibration. 


The magnitude and frequency of the steady state pitch 


“Fipure 17 depicts the actual pitch angle of the test 
craft. Figure 18 is the same data with the steady state 
average pitch angle removed. This was done on all data to 
facilitate the measurement of the pitch oscillations super- 
imposed on the steady state value of pitch angle. 
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oscillation shown in Figure 18 is approximately .08- degrees 
and .295 Hz respectively. The steady state magnitude of 

the roll motion shown in Figure 16 was approximately .075 
degrees and its frequency was approximately the ES ere Gener 
of the Deakich motion, «295 adlz. 

These steady state oscillations appear ouite small until 
it is realized that the maximum roll perturbation that was 
able to be induced during any run was only about .4 degrees. 
Lower ratios. between induced and background oscillations 
were the norm for both roll and pitch induced motion. The 
background oscillations in pitch and roll therefore could 
not be ignored when analyzing data, and deleteriously 
affected the measurement of decay time constants. 

As stated before the attempt was made to introduce 
Separately pitch and roll motion into the test craft and 
then record the transient response after the forcing action 
had passed or was stopped. After studying the processed 
data it was evident that it had been nearly impossible to 
Sancermimanements invene variahllem(either prechsor roll) 
without introducing at least some oscillations in the other 
variable. In introducing pitch oscillations roll motion was 
greatly affected. This was caused because the wave front 
caused by the chase boat's wake always struck the test craft 
at an oblique angle and not perpendicular to the axis of 
symmetry of the craft. Only during those runs where the 


XR-3 was sallied was some modicum of success achieved in 
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isolating roll transient motion. The one run included in 
this report where sallying was used to induce roll transients 
is the one conducted at 11.79 knots. During the other two 
runs included (15.93 and 19.45 knots) the wake crossing 
method Wace sed. TOm.eOrce Osetra Ons.  lhlS Meoml ted in 

runs in which both roil and pitch transients could be 

studied at the same time. 

In making comparison between the recorded data and that 
of the computer simulations both the original XR-3 model as 
presented in Ref. 3 and the new model which incorporates 
the new fan maps, new bow and stern seals and a water slope 
Correction were used. This last modification 1s One which 
has been incorporated into the XR-3 program by LCDR G. T. 
Forbes, USN during his concurrent study and is reported in 
Ref. 5. The water slope modification is used to correct 
the water level in the plenum chamber based upon the fact 
that increased pressure in the plenum will tend to displace 
water after it comes in contact with the bow seal and passes 
into the chamber. This displacement since it does not take 
place instantaneously is a function of both time and pressure 
and hence results ina "slope" in the water as it passes 
through the plenum chamber. The purpose of using both the 
original and the new models was to gain an appreciation of the 
SErCecusmetm ene MOdering changes made. 

During all the runs made at Lake San Antonio with the 


XR-3 the water conditions were fairly calm, but the steady 
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State roll and pitch oscillations were always present in 
the processed data as reported previously. Figures 19 
through 35 which are the graphs for all the runs clearly 
BAeV Coico Osc lliavion an DOU Dpiltecn and roll: Again it is 
emphasized that these oscillations precluded exact measure- 
ments of the decay time constants. 

The method used to induce rol liane’ paech transient 
oscillations in the computer models was to offset the initial 
roll and pitch angle from the known steady state angles 
when the program was executed. These offsets would then 
decay as the craft returned to the steady state conditions 
for that speed. The magnitude of these offsets were taken 
from the recorded data. The point in time where this data 
was taken is marked on each graph of the test craft's 
recorded data by a dashed vertical line. 

From known loading conditions the craft's weight during 
these runs was 6050 lbs. Its center of gravity was located 
10.08 feet forward of the transom. Average plenum pressure 
was not available as this data was not eee The plenum 
pressure was assumed to be the nominal 24.86 PSF as reported 
ma Ret. 3.. 

Teco Anents a thewcease Of the orieinal [Ref., 3] XR-3 
model caused the craft to have a negative steady state trim 
angle which is readily apparent when looking at the eraphs 
of the simulation runs of the original model (see Figures 


22 through 28 for simulation runs at 15.93 knots). 
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Although it was recognized that this negative trim angle was 
not realistic the initial conditions were not varied to get 
a positive trim angle because it was desired to test each 
model with exactly the same physical characteristics. The 
actual value of the steady state trim was not the quantity 
that was to be studied. The quantity that was of interest 
was the transient response of the model as it arrived at 


- 


that steady state value after an initial offset. 


B. RUN DATA ANALYSIS 

The first run was a calm water transient roll motion 
Buuday ao 11.79 knots. In this run the pitch offset was 
considered negligible. The roll induced oscillation was 
measured to have a negative magnitude of approximately 
.-3 degrees. Tnis roll offset was entered into both models 
With the following results. The decay time constant for 
roll motion was 4.6 seconds for the original XR-3 model and 
5.6 seconds for the new model. The oscillation frequency 
was .752 Hz for the original XR-3 and .714 Hz for the new. 
Measured data gave the actual boat a frequency of roll 
oscillation of .29 Hz and a decay time constant of less 
aati oecond. 

The second run was at 15.93 knots in calm water. This 
run had both roll and pitch perturbations. They were 
minus .25 degrees for pitch and a negative .27 degrees for 
roll. With these offsets the new XR-3 model gave a roll 


time constant of 4.1 seconds and a pitch time constant of 
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3.4 seconds. The frequencies of oscillation for the new 
model were .72 Hz for roll and .44 Hz for pitch. These 

same quantities were measured on the old model to be 

3.1 seconds and 5.0 seconds for the roll and pitch time 
constants respectively, and in order, .76 Hz and .50 Hz for 
the roll and pitch frequencies. The same measured quantities 
from the recorded data gave something less than one second 
for both the roll and pitch time constants and .25 Hz and 

.c2 Hz for the oscillation frequencies for roll and pitch. 

The last calm water data was taken at 19.45 knots. The 
measured pitch offset was a positive .12 degree and a 
positive .45 degree for the roll offset. Again these 
values were imposed upon the initial steady state trim 
conditions to achieve a tranSient reSponse. 

From the two computer simulations the following results 
were Obtained. Original model roll frequency — .76 Hz, 
roll decay time constant — 2.7 seconds, pitch frequency — 
.50 Hz, pitch time constant — 4.6 seconds. New model roll 
frequency is .72 Hz and decay time constant 4.3 seconds. 
Pitch frequency is .42 Hz and decay time constant is 
4.4 seconds. The same data from the recorded tapes gave a 
pitch frequency of .28 Hz and with decay time constant of 
less than one second. Roll frequency was measured to be 
.28 Hz, and the roll decay time constant was estimated at 


less than one second. 
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From the above data taken on the three runs and the 
computer simulations thereof it is evident that the computer 
model neither in the original form nor in its new form 
accurately predicts roll or pitch transient motion under 
the conditions studied. Decay time constants for pitch 
and roll were generally four’ times too long in both simula- 
tion models. Pitch and roll frequencies of the computer 
models were generally too high by a factor of two. Table 
IiI gives this transient behavior of the XR-3 and both 
Simulation models in tabular form. 

Figure 19 is the actual roll response data for Run 1. 
collected at 11.79 knots. Figures 20 and 21 are the 
computer simulation runs for the original and new model 
for this run depicting roll response. Recorded pitch and 
roll motion for Run 2 (15.93 knots) is given in Figures 22, 
23 and 24 and the original model computer simulation of 
roll and pitch is in Figures 25 and 26 respectively. The 
response of the new model for Run 2 is recorded in graphs 
27 and 28. Recorded data for Run 3 is given in Figures 29 
through 31 with the computer model's simulations of pitch 


and roll given in Figures 32 through 35. 


C. SIMULATION STUDIES OF THE TWO MODELS 

In an effort to gain better knowledge of the total 
effect of the new modeling, both the old and the new models 
Were given a positive three degree pitch offset from steady 


State calm water conditions. Both models were set at the 
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TABLE III. CALM WATER FORCED OSCILLATIONS RESPONSE 


nue bs opeed 11,79 knots 


Measured 
Variable XR3 
Roll Frequency .o9 Hz 
Roll decay time constant <1 sec 


Pitch Frequency - 
Pitch decay time constant ~ 


Run 2, Speed 15.93 knots 


Measured 
Variable XR3 
Roll Frequency Ay-45 te 7 
Roll decay time constant <1 sec 
Pitch Frequency eons. 
Pitch decay time constant <1 sec 


Run 3, Speed 19.45 knots 


Measured 
Variable XR3 
Roll Frequency .28 Hz 
Roll decay tine constant <1 sec 
Pitch Frequency .28 Hz 
Pitch decay time constant <1 sec 


(ith 


Simulated 


Original Model New Model 
~ foe HZ 714 Hz 
4.6 sec 5.6 sec 

Simulated 

Original Model New Model 
.76 Hz 72 Hz 
3.1 sec 4.1 sec 
50 Hz 44 Hz 
5.0 sec Bad! 

Simulated 


Original Model New Model 


.76 Hz .72 Hz 
2.7 sec 4.3 sec 
50 Hz 42 Hz 
4.6 sec 4.4 sec 
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FIGURE 19 
RoLL TRANSIENT VERSES Time, VeLocity 11,79 KNoTs 
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FIGURE 20 
Computer SIMULATION RoLL ANGLE verses Time, 11.79 knots 
ORIGINAL XR-3 
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ComPUTER SIMULATION RoLL ANGLE VERSES Time, 11,79 Knots 
New XR-3 
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FIGURE 22 
PITCH VERSES Time, VeLocity 15,93 knots 
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FIGURE 23. 
PITCH TRANSIENT VERSES Time, VecLocity 15,93 kNoTS 
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FIGURE 24 


RoLL TRANSIENT VERSES TIME, VeLocitTy 15,93 KNoTS 
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FIGURE 25 
Computer SIMULATION PitcH ANGLE VERSES Time, 15,93 kNoTS 
ORIGINAL XR-3 
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FIGURE 26. 
CompuTER SIMULATION RoLL ANGLE VERSES Time, 15,93 kNnoTS 
ORIGINAL XR-3 
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FIGURE 27 
CompuTER SIMULATION PITCH ANGLE VERSES TIME, 15,93 KNOTS 
New XR-3 
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FIGURE 28 
Computer SIMULATION RoLL ANGLE VERSES Time,15,93 KNOTS 
New XR-3 
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FIGURE 29 
ROLL TRANSIENT VERSES Time, YVetocity 19,45 knots 
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PITCH VERSES TIME, VeLociTy 19.45 KNoTs 
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FIGURE 32 
COMPUTER SIMULATION ROLL ANGLE verses Time, 19.45 knots 
ORIGINAL XR-3 
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FIGURE 33 
CompuTER SIMULATION PitcH ANGLE veRSES Time, 19.45 KNoTS 
ORIGINAL XR-3 
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FIGURE 31 
Computer SIMULATION ROLL ANGLE VERSES Time, 19.45 xKNoTs 
New XR-3 
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FIGURE 35 
Computer SIMULATION PITCH ANGLE VERSES Time, 19,45 KNoTS 
New XR-3 





same initial conditions that they were during the 19.45 knot 
run except for the offset. In the case of the new model 

this raised the bow approximately six inches and since the 
Steady state condition gave the bow a depth of about three 
inches with the positive steady state pitch angle this caused 
air gaps to appear beneath the bow seal. In the case of the 
Old model where the steady state pitch angle was negative 

no air gaps were created. Figures 36 through 38 present the 
original model's response and Figures 39 through 41 present 
the new model's response. 

As can be seen in Figure 40, the air gap on the new model 
caused wide fluctuations in the plenum pressure during the 
first half-second. This loss of pressure permitted the boat 
to rapidly sink as shown in Figure 41 until the plenum 
pressure built up. When the air gaps closed and the boat 
continued to settle the plenum pressure rapidly rose as the 
plenum volume decreased. The boat exponentially returned 
to steady state draft levels as the pitch motion shown in 
Figure 39 died down. The pitch motion of the new model 
shows the effect of the new seals. The downward movement 
of the bow as seen in Figure 39 is severely limited as the 
bow seal wetted length becomes greater and the effect of 
Bacomrieta body analysis comes into play. The first negative 
overshoot is limited to 34 percent of the original effect. 
The initial decay time constant is 2.5 seconds and the 


pitch frequency is ae Hz which because of the nonlinearity 
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FIGURE 36 


CoMPUTER SIMULATION 3 Decree Pitcu Step Response, 19,45 knots 
ORIGINAL XR-3 PITCH VERSES TIME 
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FIGURE 37 
| Computer SIMULATION 3 Decree PitcH Step Response, 19,45 knots 
Or1GINAL XR-3 PLENUM PRESSURE VERSES TIME 
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FIGURE 38 
Computer SIMULATION 3 Decree PitcH Step Response, 19.45 knoTs 
ORIGINAL XR-3 Draft at C.G. verses TIME 
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FIGURE 39 
Computer SIMULATION 3 Decree PitcH Step Response, 19.45 Knots 
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FIGURE 40 
ComPuTeR SimuLaTion 3 Decree Pitcu Step Response, 19,45 knots 
New XR-3 Prenum Pressure VERSES TIME. 
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FIGURE 41 
CompuTeR SIMULATION 3 Decree PitcH Step Response, 19.45 KNoTS 


New XR-3. Drart at C.G. VERSES TIME 
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of the bow and stern seal models changes to a frequency of 
-42 Hz at the end of ten seconds. The time constant has 
lengthened to 4.3 seconds in the same period. 

The behavior of the old model is different. Because 
of the high capacity of the fan maps the damping effect of 
plenum pressure is minimal. The draft movement shown in 
Figure 38 is therefore highly underdamped whereas the draft 
motion for the new model did not display any oscillatory 
tendencies (see Figure 41) and displayed exponential decay 
type response when returning to steady state conditions. 
In pitch motion shown in Figure 36 the original model 
displays a first overshoot magnitude of 78 percent of 
Sripiieamettser., Ihe frequency of oscillation is .5 Hz 
and the decay time constant is 4.78 seconds. Both the time 
constant and pitch oscillation frequency remain constant 
Pareueniougerme 16 Second Simulation. 

To further evaluate the differences between the original 
and new models a second simulation run was made at 19.45 
Ors. This run was in ahead sea state 2 irregular seas. 
Figures 42 and 43 show that the plenum pressure fluctuations 
in the new model were 600 percent greater than in the 
original model. Whereas the plenum pressure varied from 
5 to about 45 PSF in the new model the pressure only varied 
from about 21 to 27 PSF in the original. Draft fluctuations 
were about 300 percent greater in the new model than in the 
original (see Figures 44 and 45). Pitch angle fluctuations 


were comparable in both models (Figures 46 and 47). 
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FIGURE 42 
ComPUTER SIMULATION PLENUM PRESSURE VERSES Time, 19.45 KNOTS 


New XR-3, Sea State 2 
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FIGURE 43 
CompuTeR SIMULATION PLENUM Pressure VERSES TIME, 19.45 Knots 
ORIGINAL XR-3, Sea StaTE 2 
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FIGURE 44 
Computer SimucaTion Drart at C.G. verses Time, 19.45 knots 
New XR-3, Sea State 2 
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FIGURE 45 


Computer SIMULATION Drart at C,G, verses Time, 19.45 knots 
Or1GINAL XR-3, SEA State 2 
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FIGURE 46° 


CompuTER SIMULATION PITCH ANGLE VERSES Time, 19.45 KNoTS 
New XR-3, Sea State 2 
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FIGURE 47 


Computer SIMULATION PITCH ANGLE VERSES Time, 19.45 KNoTS_ 
ORIGINAL XR-3, SEA STATE 2 
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V. CONCLUSIONS AND RECOMMENDATIONS 


Making accurate measurements from the collected data 
was extremely difficult because of the high relative level 
of steady state oscillation present in both the pitch and 
roll motion. Signal noise coupled with the fact that the 
decay time sere Gees and the periods of the oscillations 
were of the same general magnitude made accurate measurement 
of the decay time constants practically impossible. 

Although accurate measurement of the time constants was 
not possible general observations could be made. It was 
clearly evident that the pitch and roll transients induced 
in the XR-3 test craft decayed very rapidly and not with 
the highly underdamped characteristics shown by the computer 
Simulation models. The measured tape data confirmed the 
observations made aboard the XR-3 while conducting these 
runs that the craft come to steady state conditions with 
only one or two overshoots after the forcing action had 
ceased. Roll motion transients were normally completed in 
less than two cycles. Both the longitudinal and athwartships 
Stability were highly speed dependent gaining stability as 
speed increased. 

These observations were not based solely on "feel" but 
were made by sighting across the deck of the XR-3 and 
observing the motion of the deck edge relative to the shore 


line. 
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Simulation runs have shown that for large perturbations 
the transient response of the new model is quite different 
from the original. Since it was not possible to excite such 
large perturbations in the actual test craft it remains to 
be verified that the new bow and stern seal models more 
accurately predict the forces generated by the actual seals 
than does the original seal model. 

The new seal fan maps are more accurate than the original 
Since they are based on actual verified manufacturer's data. 
' Leakage areas are now proportional to those of the actual 
Clean . 

The fact that the present computer models of the XR-3 
both old and new do not appear to accurately represent the 
actual transient behavior of the XR-3 in pitch and roll 
opens many areas which should be investigated for possible 
errors. AS both the pitch and roll transients appear to be 
in error the sidewall subroutine becomes a prime area of 
investigation as it contributes major components to both 
the roll and pitch force moments. 

Figure 48 shows one area which should be explored and 
that is the area of hydrodynamic lift of the sidewalls. 
Figure 48 is a simplified drawing of the underside of 
the XR-3. The flat sections on the bottom of each 
sidewall are cross-hatched for clarity. At high 
velocity (20-25 knots) these areas become effective 
planing surfaces. In this program the lift of the 


Sidewall sections due to planing does not appear to take 
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into account the underwater body shape. This shape changes 
radically from bow to stern. Whether or not this shape has 
Significant effect on roll and pitch transient behavior is 
not known but should be investigated. Another aspect of 
the sidewall modeling which should be investigated is the 
inclusion of forces due to the displacement of water per 
unit of time. A term such as this would not only have terms 
based on draft and pitch angle but would also include terms 
which took into account the rate of change of displacement 
and the rate of change of the pitch angle. The addition of 
first derivative terms such as these would have a definite 
effect on the transient behavior of the sidewalls. 

A third area of study concerning the sidewalls is the 
effect of sidewall deadrise on turning performance. Modeling 
of the deadrise forces would be incorporated in the cross- 
flow force modeling. This deadrise could have significant 
effect on the outboard sidewall during a turn determining 
whether it has a tendency to rise and ride over the water 
as the craft skids or dig in and cause the craft to undergo 
Violent roll motion as reporeed ain Ref. 6. 

Another area of study which needs investigation is the 
effect of plenum pressure on the transient response (of the 
computer model). If it can be shown that the computer 
model as presently configured can be made to duplicate 
actual measured craft transient motion by just adjusting 
View pienum pressure, Turther modifications of the program 


may be unnecessary. Further transient motion studies are 
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needed where both the plenum pressure and the seal pressures 
are included in the recorded data. Because of the stochastic 
nature of the actual measured data many duplicate runs will 
be needed to arrive at accurate averages of the test craft's 
oscillation frequencies and the decay time constants. These 
averages should then be used in comparing the actual craft 
with the computer model. 

As ‘canbe seen in@the recerded data presenved herein a 
fruitful area of study would be in reducing the signal noise 
on the recorded data. This could be in the recording phase 
or processing phase of the data. Thorough tests under 
controlled conditions need to be made of the actual transfer 
Characterisuics of the different Sensors installed in the 
wocmecralt, 

Earlier in this report it was mentioned that there 
existed steady state roll and pitch oscillations while the 
boat was traveling straight on calm water. These oscillations 
could exist and could be the result of the craft being in 
some limit cycle type of motion because of the inherent 
nonlinearity of vehicles operating in fluid mediums. On 
Vaemornertmnand. it could be Ehesresult sor instrumentation 
for example if the vertical seeking gyro from which all roll 
and pitch motions are derived was precessing about the true 
vertical then it could impose a sinusoidal oscillation 
signal on both the pitch and roll sensors that would show 


up in the recorded data. 





In conclusion it must be stated that there are still 
many areas in this subject of sidewall modeling of the 
XR-3 which should and need to be explored before the XR-3 
model as originally stated in Ref. 3 and modified herein 
can be accepted as a reasonably accurate representation of 
the actual XR-3 test craft. Further investigation into 
the data collection and processing is also needed to refine 
the processed data to the point where it is accurate 


enough for dynamic studies. 
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APPENDIX A 
LISTING OF BOWSEAL SUBROUTINE ‘BOWSL 


The entering arguments needed for the new bow seal 
routine are the same as presented in Ref. 3 for the original 
bow seal except that instead of entering the angle between 
the leading edge and craft vertical in colums 46-55, the 
length of the rear cables is entered (in feet). An additional 
card is needed to enter the length of the middle support 


cables. The new data card should be entered as shown below. 


CARD COLUMNS FORMAT ENTRY 
il ~}- 5 I 00601 — control tag 
6-15 F X coordinate of seal 


hinge (ft) forward 
of transom 


16-25 F Seal leakage orifice 
coefficient 
26=35 F Pressure differential 


between bow seal and 
plenum chamber (PSF) 


36-45 F Z coorginate of seal 
hinge (ft) above keel 

46-55 Length of rear support 
cable (ft) 
(17675 “fs full down) 

56-65 Length of the seal 
leading edge (ft) 

66-75 Base leakage area (ft*) 
Nominal .08 ft°) 

2 1- 5 00602 — control tag 
6-15 Length of middle support 

cable (ft) 


(1.1875 is full down) 


IES. 
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APPENDIX B 
LISTING OF STERN SEAL SUBROUTINE ‘STNSL 


The entries for this routine are the same as the original 
stern seal except that in columns 46-55 the rear support 
cable length is entered instead of the seal angle. 


Block 5 entries are as follows: 


CARD COLUMNS FORMAT ENTRY 
il: l= 5 Tt 00500 — control tag 
6-15 F X coordinate of seal hinge, 
(ft) forward of transom 
16-25 F Z coordinate of seal hinge 
above keel (ft) 
26-35 F Base leakage area (f°) 
See Omei a 
36-45 F Seal leakage orifice 
coefficient 
46-55 F Length of rear support 
Gomme (ft) 1.875 is full 
down 
56-65 F Pressure differential 


between stern seal and 
plenum chamber (PSF) 


66-75 F Length of the leading edge 
of the seal (ft) 
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APPENDIX C 
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LISTING OF FAN SUBROUTINE FAN 


The entries for the fan subroutine remain the same. 


Printed 


The nominal fan RPM is 8000 RPM vice 2615 RPM. 


below are copies of the appropriate entries for the bow 


stern seal and main plenum fans. 
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APPENDIX D 


LISTING OF INITIAL CONDITION SUBROUTINE INCON’ 


The INCON subroutine was modified to permit it to 
calculate the wave components of the various sea states if 
desired instead of having to enter each component individually. 

Block 11 options 1 and 2 remain identical. New options 
3 and 4 have been added and are entered as follows: 


BLOCK 11 - OPTION TAG 3 


CARD COLUMNS FORMAT ENTRY 
AL 1- 5 I OTIOS Control and option vac 
6-15 F Number of wave components 


maximum of 10 


16-25 F Initial wave heading (deg) 
(180 deg = ahead waves) 

26-35 F Average height of 1/3 
highest waves (ft) 

36-45 F Shortest significant period 
(sec) 

H6-55 F Longest significant period 
(sec) 


BLOCK 11 - OPTION TAG 4 


CARD COLUMNS FORMAT lau Alaé 
1 1- 5 I Olli weontrol and Option tag 
6-15 F Number of wave components 


maximum of 10 


16-25 oor Initial wave heading (deg) 
(180 deg = ahead waves) 


26-35 F Average height of 1/3 
highest waves (ft) 


ees, 





CARD COLUMNS FORMAT ENTRY 


1 36-45 F Lowest significant frequency 
(Hz) 

46-55 F “ae significant frequency 
Hz 
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